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I— I Abstract: We study dijet production in proton-proton collisions with a veto on the 

p emission of a third jet in the rapidity region in between the two leading ones. We 

p resum the leading logarithms in the ratio of the transverse momentum of the lead- 

^ ing jets and the veto scale and we match this result to leading-order QCD matrix 

I— I elements. We find that, in order to obtain sensible results, we have to modify the 

y—^ resummation and take into account energy-momentum conservation effects. We com- 

^ pare our theoretical predictions for the gap fraction to experimental data measured 

OO by the ATLAS collaboration and find good agreement, although our results are af- 

i^ fected by large theoretical uncertainties. We then discuss differences and similarities 

(^ of our calculation to other theoretical approaches. 

O 
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1. Introduction 

QCD phenomenology and in particular jet physics are playing a central role in the 
physics programme of the first period of the LHC running [1-7]. While the predictions 
of inclusive jet and dijet cross sections are under good perturbative control, the 
prediction of the associated final states is much more delicate. A crucial role is played 
by the colour structure of the hard process, which sets the initial conditions both of 
the parton shower, which describes the perturbative evolution of the event and the 
hadronization, which describes the transition of those partons into the hadrons that 
are seen in the final state (for an overview of this physics, see for example [8]). While 
long viewed as a probe of QCD dynamics, it has only more recently become widely 
realized that a measurement of hard process colour structure could be an important 
probe of new physics [9-12] (although the idea dates back to at least [13]). 

The most direct probe of the colour structure of a hard process is the probability 
that it does not radiate into some well-defined region of phase space. In this paper 
we consider the prototypical process, in which a dijet system is measured and the 
presence of any additional jets in the rapidity interval between them is vetoed. Pro- 
vided the threshold for reconstructing and vetoing these jets is in the perturbative 
regime, the cross section for such events can be calculated perturbatively. However, 
at every order of perturbation theory, logarithms of the ratio of the hard process 
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scale to this jet veto threshold scale arise, and when this ratio is large it is manda- 
tory to sum these logarithms to all orders. When the ratio of scales is not large, these 
logarithms do not dominate and fixed order perturbative results are more accurate 
than the resummed results. To provide a prediction that is valid over all values of 
the ratio, it is necessary to match the resummed and fixed-order calculations, which 
is one of the aims of this paper. The jet veto "gap fraction" has been measured by 
the Tevatron and HERA collaborations [14-18] and, very recently, for the first time 
at the LHC [3]. 

In a previous paper [19] two of us made a first phenomenological study of this 
observable, comparing the resummation of soft gluons to a standard event generator 
(Herwig-|--|- [20]). We found that neither approach was completely satisfactory. The 
parton shower simulation does not contain the full colour structure or any account 
of the mixing between different colour structures, and neglects contributions coming 
from loop corrections that do not correspond to a non-emission contribution (so- 
called Coulomb gluon exchanges). These give a sizeable correction at large enough 
jet transverse momenta and rapidity separation. Moreover, higher than expected 
non-emission probabilities were found, which has more recently [21] been explained 
as a problem with the way Herwig-|--|- implements the colour structure of hard 
processes involving gluons. On the other hand, the resummation approach, which is 
based on the soft gluon approximation, by itself is not sufficient for currently-probed 
phase space regions. One important effect that it neglects is energy-momentum 
conservation: within the soft gluon approximation, the additional jets that one is 
vetoing could be emitted with a suppression only due to matrix elements, whereas in 
full QCD such emission would be considerably suppressed by the requirement that 
it carry away enough energy that it be above the jet veto scale. 

In this paper we improve the previous resummed predictions for the gap fraction 
by modifying them to approximately account for energy-conservation effects and by 
matching them to the leading QCD order calculation. We also include the first tower 
of non-global logarithms arising from one gluon emission outside of the gap region. 

As well as providing an important understanding of the QCD effects that de- 
termine emission patterns, a precise calculational framework would enable colour 
structures to be used to constrain the physics processes leading to new particle pro- 
duction. This has been particularly thoroughly studied for the Higgs boson in asso- 
ciation with two jets [12,22,23], for which the colour structure has been shown to be 
identical to the dijet processes that we study here [24]. 

In the remainder of this paper, we define more precisely the observable we will 
be calculating, the gap fraction. We then discuss the all-orders resummation of 
the associated leading logarithms, and how to match this result with a fixed-order 
result. Finally, we compare the matched results with the ATLAS data and with 
other theoretical predictions, before drawing some conclusions. 
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2. The gap fraction 

We are interested in dijet production in proton-proton collisions: 

where we veto on the emission of a third jet with transverse momentum bigger than 
Qq in the rapidity region between the two jets. In the present study we fix the 
veto scale at Qo = 20 GeV. Pi^2 define the incoming hadron momenta and ^3^4 the 
outgoing jet momenta. We define the gap fraction as the ratio of the cross section 
for this process over the inclusive rate: 

^2 gap , J2 
/gap ^ ^±_Z__ /____ (2 I) 

■' dQdY/ dQdY' ^ ■ ^ 

In the Born approximation, the final state consists only of the two hard jets, so 
every event is a gap event and /^^^ = 1 . Beyond the Born approximation the leading 
jets are no longer balanced in transverse momentum. We define Q to be the mean 
of the transverse momenta of the leading jets Q = (prs + Pta)/'^- This choice, in 
contrast for instance to the transverse momentum of the leading jet, is more stable 
under the inclusion of radiative corrections [25]. The rapidity separation is defined 
by y = Ay — 2D, where Ay = {y^ — y^l is the rapidity separation between the centres 
of the leading jets and D can be freely chosen. In many previous studies, D was set 
equal to R, the jet radius. The ATLAS collaboration instead measure the gap region 
from the centres of the jets, i.e. D = and thus Y = Ay. It is useful to rewrite 
Eq. (2.1) using unitarity: 

dVs^P dV^^ _ dV 

dgdF^dgdF ~ dgdy' ^^'^^ 

where we have introduced the complement of the gap cross section, which corresponds 
to requiring at least one jet harder than Qq in the rapidity region in between the 
leading ones. The gap fraction then becomes 

dV^^ / dV 
^'^'^^~ dQdY/ dQdY- ^^-^^ 

Our target in this paper is to evaluate this expression at the first order in the strong 
coupling (LO) and to match it to the resummed calculation. At this accuracy then 
Eq. (2.3) contains only tree level contributions; the numerator is the integrated trans- 
verse momentum distribution of the third jet over the gap region and the denominator 
is simply the Born cross section: 

'" ^ ' - L *^ I * ^^ d..d,d^dQdy / dQdy + °("2'- (^■^' 
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where kx, y and are the transverse momentum, the rapidity and the azimuth of the 
third jet. The notation J. imphes the integral over the gap region in rapidity and 
azimuth. We compute the gap fraction for proton-proton colhsions at 7 TeV using 
NLOJET++ [26] at leading order. The jets are defined using the anti-kt algorithm [27] 
with R = 0.6. We use the Cteq6.6 parton distribution functions (PDFs) [28] and 
adopt the same kinematical cuts as the ATLAS collaboration, requiring all jets to 
have pt > 20 GeV and \y\ < 4.4. 

A fixed-order calculation of the gap fraction is reliable only at small Ay and 
when Q is of the same order as Qo- As soon as we move away from this region, 
the leading-order gap fraction decreases rapidly and eventually becomes negative. 
This unphysical behaviour indicates that the fixed order calculation by itself is not 
reliable. Large logarithms of the ratio Q/Qq contaminate the perturbative expansion 
and they must be resummed to all orders, as discussed in the next section. Also 
terms proportional to Ay (formally equivalent to a logarithm) can be resummed, for 
instance as in the High Energy Jets (HE J) framework [29] . In the limit of large Ay 
and Q/Qo, the cross section is dominated by the singlet exchange component and 
there is overlap between the logarithms resummed by the approach we describe here 
and those resummed by the BFKL equation [30]. 



3. The resummed calculation 

The technique for resumming logarithms of the ratio Q/Qo for the gaps-between-jets 
cross section has been explained in detail in [19,31-35]. It relies on the ability to 
map the real part of loop corrections into a form that is exactly equal and opposite 
to the phase space integral for real emission. For an observable in which emission is 
suppressed equally in all angular regions of the event, a global observable, there is 
an exact cancellation between the real and virtual contributions, such that the result 
(up to a phase term) is a virtual integral over the part of phase space in which real 
emission is vetoed. 

In the observable we are studying, however, radiation is only suppressed in part 
of the phase space region and not globally. For this reason "in-gap" virtual correc- 
tions are not enough to capture even the leading logarithmic accuracy. Radiation 
outside the gap is prevented from re-emitting back into the gap by the veto require- 
ment, inducing a further real-virtual miscancellation and additional towers of leading 
logarithms, called non-global logarithms [36]. Currently, these contributions can be 
resummed only in the large Nc approximation [37,38]. Here, instead we adopt the 
approach suggested by [34]: we keep the full colour structure but we expand in the 
number of gluons, real or virtual, outside the gap. It was argued in [19] that this 
may be a reasonably convergent expansion, so that the full result is approximated 
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by the contributions arising from only zero or one gluons outside the gap: 

^2 gap ^2 (0) ^2 (1) 

dQdY dQdY dQdY ^ ' ' 

The first contribution to this expansion corresponds to the exponentiation of the 
one-loop virtual corrections (with no gluon outside the gap). 
We define the resummed gap fraction as 

J2 gap ,(^2 born 
^gap ^ ^ '^rcs / ^ " (o 2) 

■'''' dQdY/ dQdY' ^ ■ ^ 

Because we are working in the eikonal approximation, additional radiation does not 
change the Born kinematics and the resummed cross section factorizes into products 
of resummed partonic contributions and parton luminosity functions: 



d^cT^*) P^raf 



E TX^TX^I-^2c.r>C„,(A|/,g) (3.3) 

^-^ . 1 + dnh 1 + Ord A?y=Y 



dQdr 2QS ^l + 6abl + Scd 

a,b,c,d 



with 



CatiAy, Q) = ^ f ^ dy fa{^zey'\ Q)h{^ze-y'\ Q) , (3.4) 



where z = xiX2, where xi 2 are the longitudinal momentum fractions of the incoming 
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partons, y = 1/3 + Vi and p = ^. The integration limits are 



I/+ =min ( ln-,|/cut j , (3.5) 

where the value ijcnt is obtained by requiring that both jets are within the calorimeter 
acceptance. 

The resummation of global logarithms is achieved by considering the original 
four-parton matrix elements dressed by in-gap virtual gluons, with transverse mo- 
menta above Qq and no out-of-gap (real or virtual) gluons. The resummed partonic 
cross section then has the form 

|A^(°f = tr (^ife-^^Q^'^^^'^e-^^Q'-Q^r) , (3.6) 

where ^ is computed by considering the strong coupling at one loop: 

,(1. l.^ 2 f'^ dkr ,, , 1 l + «.(g)/3oln| 

^{ki,k2) = - -—as{kT) = —^\n — ^, (3.7) 

with (3o = — ^^^ "^ . The matrix H in Eq. (3.6) gives the matrix elements of the hard 
process in some colour basis, while S is the metric tensor in that colour basis. In an 
ortho normal basis, as we use throughout this paper, S = 1. 
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The soft anomalous dimension is a matrix, with elements Fjj = (ej|r|ej) in a 
basis {cj}: 



2 L 2/'^ 



111 



ta ■ tf, Co'12 + ta ■ tc W13 + t^ ■ trf ^14 



+tfe ■ tc 6^23 + tf, ■ td W24 + tc ■ t(i iOs4 



+ irrta ■ t 



■ ^b 



where tj is the colour charge of parton i^. The soft gluon momentum is labelled 
k. The integrals over the gluon's azimuth and rapidity inside the gap region admit 
simple analytical expressions if one considers azimuthally symmetric gaps. However, 
in the current analysis, we are defining the gap from the centres of the leading jets. 
Thus, the gap region is not just a rectangle in the [r], 0) plane and we have to integrate 
around the two semi-circular boundaries of the leading jets. Analytical expressions 
can still be obtained as a power series in R [39] . Here instead we decide to keep the 
full R dependence and perform the integrals numerically when we cannot find simple 
analytical results. For the explicit expressions of the hard scattering matrices H in 
the various partonic channels we refer to [19]. 

3.1 Non-global contribution 

We want to estimate the impact of non-global logarithms on the gap fraction. In par- 
ticular we aim to resum the non-global logarithms that arise as a result of allowing 
one soft gluon outside the rapidity gap. The general framework in which this calcula- 
tion is performed is described in [34,35], where the case of an azimuthally symmetric 
gap was considered. As in the global case this led to relatively simple analytical 
expressions. To include the gap definition used in the ATLAS analysis instead, one 
has to resort to evaluating most of the integrals numerically. This considerably slows 
down the calculation, but has a very small effect on the final results. Therefore, for 
the current work, we decide to include non-global effects as a K-iactoi: 

T^/r^ A \ dQdY ^ dQdY /„ n^ 

K{Q,Ay) = ^,^(0) , (3.9) 

dQdF 

where in calculating this ratio we compute the resummed cross sections for az- 
imuthally symmetric gaps, and then use it to multiply the resummed result for zero 
gluons outside the gap including the exact gap definition. The error this approxima- 
tion induces is much smaller than the overall uncertainty in the resummed approach. 



^Note that the mismatch between the indices 1,2,... and a,b, . . . is related to the fact that one 
must sum over two orientations of the event, e.g. i{pi,ta) + j{p2,'tb) -^ ^(P3,tc) + l{p4,td) and 
i{pi,ta) + j(p2,tb) -^ l{p3,td) + fc(p4,tc), as explained in more detail in [19]. 
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which we estimate below. This approximation does not affect our matching proce- 
dure because we are only performing LO matching and non-global logarithms start at 
O (a^) in the expansion of the gap fraction. We do not need to include effects related 
to parton recombination due to the particular choice of the jet-algorithm [40,41] since 
we employ the anti-kt jet algorithm [27]. 

The calculation of the contribution from one gluon outside the gap is essentially 
that presented in [34,35], except that the final integral over that gluon's momentum 
is explicitly performed numerically: We briefly recap the results. The soft anoma- 
lous dimension for four-parton evolution in the case of azimuthally symmetric gaps 
reduces to 

r = -re, + i-nK ■ t, + -p{Y- \^y\){tl + t^), (3.10) 

where 

, , , , sinh (Ay/2 + W2) 

and tt = ta + tc is the colour charge matrix corresponding to emission from the total 
colour exchanged in the t channel. To obtain the contribution from one gluon outside 
the gap we must now consider both real and virtual corrections to the four-parton 
scattering, each dressed with any number of soft gluons: 

|-M('f = — / T^«,(M / di/(fiK + f]y), (3.12) 

where the integrals are over the transverse momentum and rapidity of the real or 
virtual out-of-gap gluon. The operator to insert this gluon off the external legs is 

D'' = taK + tbhi^ + t,h^ + tahl h>t=\kT^, (3.13) 



2 Pi-k- 



for real emission, and 
Ir 



1 = — 
' 2 



ta ■ tfe UJi2 + tfl ■ tc Wi3 + ta ■ td Co'14 + tf, ■ t^ UI23 -|- tfe ■ trf aJ24 + tc " trf U34 



Uij = -4 ^]'P' , , (3.14) 

for virtual emission. In the case of the out-of-gap gluon being virtual, the subsequent 
evolution is unchanged from that of the original four-parton system, given by T in 
Eq. (3.10). In the case of real emission we have to consider the colour evolution of a 
five-parton system [42,43]. If we assume the gluon to be emitted on the same side 
of the event as partons a and c, the anomalous dimension is given by 



A = ^FT? + ^vrT, ■ T, + \p{Y- Ay){Tl + T^) + lp(F; 2 y\)Tl 




+ \kY- Ay , |2/U)Tc ■ Tfc, 


(3.15) 


Tl = {T, + Tdf 


(3.16) 
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Figure 1: The i^- factor we use to estimated non-global effects, as a function of Q for 
different values of Ay (on the left) and as a function of Ay, for different values of Q (on 
the right). 



and we have introduced the kinematic function 



1 cosh( Ay/2 + y + Y)~ sgn{y) cos 

KY; Ay, y,(t>) = t; In — , /. /^ , v~\ r~\ J " ^■ 

I cosh(Ay/2 + y — y ) — sgn(|/j cos 



(3.17) 



The real and virtual out-of-gap emissions, dressed to all orders with in-gap virtual 
corrections are thus 



^R = tr 
fly = tr 



'^g-5(fcT,Q)rtj-)^tg-5(Qo,fcT)Atg-g(Qo,fcT)Aj3 g-5(fcT,Q)r 



'^g-5(Qo,Q)rtg-5(Qo,fcT)r^g-5(fcT,Q)r _^ ^_^_ 



(3.18) 



The i^-factor defined in Eq. (3.9) is plotted in Fig. 1 as a function of Q for 
different rapidity values and as a function of Ay for different values of Q. In the 
case of the Q distribution, we see that the effect is modest for the first rapidity bin, 
but is typically of the order of 30% in much of the Q range we study. Non-global 
effects are much smaller in the case of the Ay distributions for the regions of Q we 
are considering. 

It has been shown [34, 35] that naive QCD coherence is violated at sufficiently 
high perturbative orders because the Coulomb gluon exchange terms included in 
r and A induce a mis-cancelation between the real and virtual contributions in 
Eq. (3.12). The y — > oo region therefore gives a finite contribution and, as a con- 
sequence, super-leading logarithms {a^log'^~^^{Q/Qo)) arise at C(as) and beyond. 
The numerical impact of these contributions has been studied in [19] and found to 
be generally modest. 
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4. Matching 

The resummed calculations are based on the eikonal approximation, in which energy- 
momentum is not conserved. There is no recoil of the hard lines against the emissions 
and no account of their effect on the incoming parton momentum fractions at which 
the parton distribution functions are evaluated. Matching to the full 2—7-3 matrix 
elements takes into account energy-momentum conservation, at least for the first 
(hardest) emission. Its energy is taken into account and it is hence is less likely, so 
the matched gap fraction will be bigger than the one given in Eq. (3.2). 

In this section we discuss the matching of the resummed gap fraction to the LO 
calculation. At LO we can write Eq. (2.4) as follows: 



fgap 
J LO 



2as{Q) 



TT 



ao{Q,Y)\n^-bo{Q,Qo,Y) 
vo 



+ Oial 



(4.1) 



where the contribution bo is now free of large logarithms of Q/Qo- We want to 
combine this expression with Eq. (3.2), subtracting the double counted term 1 — 
^^aoln((5/Qo)- Firstly, we have to verify that the two calculations agree in the 
asymptotic limit ln(Q/Qo) ^ 1- This can be easily achieved considering the loga- 
rithmic derivative of the gap fraction: 



2as{Q) 



7T 



ao{Q,Y) = lim 



Qo^odlnQo-'^^ 



dSf^gap 



2^born 



'dV 



lim r_l /iLl f4 21 

Qo^odlngodgdF/ dQ dV ^ ' ' 



The result is shown in Fig. 2: the logarithmic derivative of the gap fraction is plotted 
as a function of ln{Qo/Q) for fixed kinematics. In this particular example we have 
Ay = 3 and Q = 200 GeV. The plot shows that the logarithmic derivative of the 
gap fraction tends to a constant for large, and negative, values of the logarithm. The 
numerical value is in agreement with the one obtained by expanding the resummation 
at 0{as{Q)): 



fgap 
J res 



l-ao{Q,Y)^-a,{Q,Y)e 

Q 



l-'-^^aMY) 



TX 



Qo ^-^. .In., kr kn 



n=l 



+o{e 



(4.3) 



The plot in Fig. 2 shows that we have control of the logarithms at 0{as). How- 
ever, plotting instead the Q dependence at fixed Qq for various Ay bins, as in Fig. 3, 
the picture is not so clear. Because the plot is on a logarithmic x axis, we might 
naively expect the FO result (data points) to asymptotically tend to a straight line, 
with the same slope as the expansion of the resummation (dashed curve). However, 
changing the Q values, one changes the momentum fractions and factorization scales 
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A y = 3, Q = 200 GeV 



FO 
Expansion 



1.5 - 
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0.5 



llHlI;JI^HlTilSTT^I^I,^^I^,I^,I, 



I /^i^K^^ 



-5 -4 -3 -2 -1 

ln(Qo/Q) 



Figure 2: The logarithmic derivative of the gap fraction ^j^^ as a function of InQo/Q 
for fixed Q = 200 GeV and Ay = 3 ("vAS = 7 TeV). The solid line is the coefficient obtained 
by expanding Eq. (3.2) at 0{as)- 



of the parton distribution functions and hence the mix of different flavour processes, 
so one could expect some curvature, but this effect should also be included in the 
expansion of the resummed results, where some curvature is also seen, so the dif- 
ferences in slope between the data points and dashed curves is really significant. 
Because the FO curve and the expansion of the eikonal resummation differ so much, 
a simple matching procedure in which we add together the FO and the resummation 
and subtract their common term, is bound to fail. It is clear that this issue must 
be investigated in more detail. The strengthening of the curvature at the highest Q 
and Ay values indicates that we are becoming sensitive to the kinematic limit and 
we therefore examine the issue of energy-momentum conservation. 

4.1 Energy- momentum conservation 

The resummed cross section (Eq. (3.3)) has been obtained in the eikonal limit, 
i.e. emitted gluons are considered soft and they do not change the Born kinematics. 
Even if we are guaranteed that this assumption is sufficient to capture the leading 
logarithmic behaviour, we are losing important physical effects related to energy- 
momentum conservation. In particular, because of the choice Qo = 20 GeV we are 
sensitive to emissions of gluons with non- negligible transverse momentum with re- 
spect to Q. Furthermore, the emission of a gluon requires a finite amount of energy 
and, for given Q and Ay, this means we are probing the parton distribution functions 



-10 



1.5 



-1.5 



Gap fraction: FO and Expansions 
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Figure 3: The gap fraction at 0{as) as a function of Q in different Ay bins. The points 
are the exact FO calculations, the dashed ones the expansions of the eikonal resumma- 
tion (Eq. (4.3)) and the solid curves correspond to the 0{as) expansion of the modified 
resummation. 



at larger values of Xi 2- Because the PDFs are steeply falling functions of x at large x 
this can give a considerable suppression well before we reach the edge of phase space. 
We would therefore like to go beyond the soft approximation and modify our 
resummation so that we can capture the correct kinematic behaviour, at least for 
the hardest (i.e. highest kx) gluon emission. In order to do that we study the 
full kinematics of a 2 — )■ 3 process and using energy- momentum conservation we 
determine the values of Xi 2: 



2^1,2 = A±e^^ , with A± 



Q cosh — ± Q sinh 1 e^^ 



(4.4) 



The variables Q, Ay and y are the same as the ones previously defined. We see that 
Xi,2 depend on the transverse momentum of the emitted gluon kr and its rapidity in 
the partonic centre of mass frame, y'. We have also introduced 



Q = Pt3- Pta 



kx kx + 2Q cos < 



2 2Q + kx cos ( 
As a consequence the parton luminosities become 



(4.5) 



Cab{Ay,Q,k) 



2A^A^ 



dyfa{A+ey/',Q)MA 



-5/2 



Q) 



(4.6) 
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with 



y~ = max {inAi, -ycnt) 
1 

A 



y+ = mill ( In ^^, t/cut ) • (4.7) 



If we take the hmit kx ^ then Eq. (4.6) reduces to the parton luminosity computed 
in the soft hmit, Eq. (3.4). 

So far we have discussed how to take into account the complete kinematics in 
the PDFs for the hardest emission. Clearly, the matrix elements will also differ from 
their eikonal approximations and the gap fraction (/ = 1 — dcr^'^P/dcr'^°™) involves 



do-g'^P 



jdkTdy'd<l)\M2^3iAy,Q,k)\'CiAy,Q,k) , (4.8) 



dgdy 

where we have suppressed parton indices for clarity. Since both the matrix elements 
and the parton luminosities depend on the momentum of the emitted parton, k, we 
lose the convenient kinematic factorization of Eq. (3.3). Importantly, it is the shift 
in the argument of the PDFs that dominates, and so we shall evaluate the matrix 
elements in the eikonal limit. To further simplify matters we can also restore the 
kinematic factorization by approximating the integral of the parton luminosity by 
its value at a particular phase space point. Specifically, we write 

drr^^ /■ 

— — ^ dkTdy'd^\M2^3{Ay,Q,k)\l,,C{Ay,Q,k) (4.9) 

dQdY J kT=VQM,y'=a^y 

and the value of Q is determined by its azimuthal average: 

. gO^-l-Ol^). (4.10) 

For the transverse momentum, we have chosen the geometric mean of the integra- 
tion limits, i.e. kr = VQoQ- The rapidity value is determined by requiring the 
approximate result on the right hand side of (4.9) to be as close as possible to its 
exact value on the left hand side. We keep a fixed as we vary Q, but allow it to vary 
with Ay and typically find | ^ a ^ | . 

The 0{as) modified gap fraction is plotted in Fig. 3 with solid lines. The plot 
shows that with a one parameter fit we can construct a modified resummation whose 
first-order expansion reproduces the FO result very accurately. We stress that this 
modification of the parton luminosity does not affect the formal leading logarithmic 
accuracy of our calculation. Rather it corresponds to a particular choice of important 
sub-leading terms that is motivated by energy-momentum conservation: it is very 
reassuring that such a procedure reproduces so well the exact leading order result. 

Using this modified resummed result we are now ready to complete the match- 
ing to LO (the matching corrections are now very small), estimate the theoretical 
uncertainty and then compare to data. 
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5. Matched results and comparison to data 

We define our modified resummed cross section (for zero gluons outside the gap), i.e. 
the replacement of Eq. (3.3), as 



a,b,c,d 



xCai>iAy,Q,k) . (5.1) 

We then define a resummed gap fraction by adding the FO calculation and the mod- 
ified resummation together, subtracting the expansion of the resummed expression 
to 0(a,): 

fgap _ ^gap I /.gap _ /.gap I'c; 9^ 

^matched ~ J LO ^ J mod Jraod,as- W'^J 

We also estimate the effects of non-global logarithms by multiplying the above ex- 
pression by the K-ia.ctoT defined in Eq. (3.9). The calculation we have performed 
matches together a LO computation with a leading logarithmic one and so we expect 
it to have a considerable theoretical uncertainty. Because we are considering the gap 
fraction, renormalisation and factorization scale variations do not give the dominant 
contribution to the uncertainty. Parton-distribution-function effects also largely can- 
cel in the ratio. The dominant source of uncertainty comes from higher logarithmic 
orders in the resummation. In particular, a leading logarithmic resummation does 
not fix the argument of the logarithms we are resumming. As an estimate of our 
theoretical uncertainty we then rescale the argument of the function ^, Eq. (3.7): 

aQo,Q)^aiQo,Q), (5.3) 

with 7 allowed to vary in a range of order unity. Motivated by the fact that next- 
to-leading logarithmic corrections to leading soft logarithms are typically found to 
be negative, we consider variations in the upward direction by, quite arbitrarily, a 
factor of 2, i.e. 1 < 7 < 2. 

The plots in Fig. 4 show the gap fraction as a function of the mean transverse 
momentum of the two leading jets Q in four different rapidity bins, while the ones 
in Fig. 5 are for the Ay distribution, in two different Q bins. We use the same 
cuts as the ATLAS analysis: all jets must have pt > 20 GeV, \y\ < 4.4 and the 
mean transverse momentum of the two highest transverse momentum jets must be 
Q > 50 GeV. The dash-dotted red line represents the LO calculation, the dashed 
green line the resummed gap fraction in the eikonal limit, solid blue is the resummed 
and matched result Eq. (5.2), with the band obtained by varying 7, as explained 
above, and finally the magenta band corresponds to the resummed and matched gap 
fraction with the non-global effects included. The black crosses are the data points 
measured by the ATLAS collaboration [3] with the gap defined by the two highest 
Pt jets (we have combined the statistical and systematic errors in quadrature). 
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Figure 4: The matched gap fraction as a function of the transverse momentum Q in 
different rapidity bins. 



The FO calculation is clearly only sensible in the first rapidity bin and for Ay > 2 
it decreases very rapidly as a function of Q and eventually becomes negative. This 
unphysical behaviour is driven by a large logarithmic term ~ as Ay In ^ which needs 
to be resummed. The eikonal resummation restores the physical behaviour but, 
as we have previously discussed, completely ignores the issue of energy-momentum 
conservation and produces too small a gap fraction. Our matched curves, with the 
inclusion of non-global logarithms, does seem to capture most of the salient physics. 
However, our results are affected by large theoretical uncertainties due to the fact 
the calculation is accurate only at the leading logarithmic level. The extension of 
resummation for the gap cross section at the next-to-leading logarithmic accuracy 
will certainly reduce this uncertainty but it is not an easy task and it is not likely 
that it is going to be completed soon. Another way of reducing the uncertainty is to 
perform the matching at NLO, so that any dependency on the rescaling factor 7 is 
pushed one order higher in the strong coupling. With the necessary NLO calculations 
available in Nlojet++ [26], such a NLO matching is certainly feasible. 
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Figure 5: The matched gap fraction as a function of rapidity separation Ay in two different 
transverse momentum bins. 



6. Comparison to other approaches and conclusions 

In Ref. [3] comparisons are made between the data and the predictions of some of 
the different theoretical tools currently available. Firstly we notice that gap frac- 
tions are defined there with respect to the dijet cross section at NLO, while we use 
the Born cross section. We have checked that, because of the definition of Q as 
the mean transverse momentum, NLO corrections are small. The best description 
of the data was found using POWHEG [44-46], interfaced with PYTHIA [47]. The 
results obtained using POWHEG interfaced with HERWIG++ [20] were found to un- 
dershoot the data. The difference between the two parton showers can be taken as 
indicative of the theoretical uncertainty due to the parton shower-NLO matching. 
The formal accuracy of the POWHEG calculations appearing in the ATLAS paper 
is not different to ours: Tree-level matrix elements are used and then matched to 
a parton shower, which is essentially a leading logarithmic resummation. However, 
the final predictions differ from ours because of the assumptions and approximations 
contained in the showering algorithm. Firstly, energy-momentum conservation is 
properly accounted for in every emission in a parton shower, not just the hardest as 
in our calculation. Also, the parton shower is limited to the large N^ approximation, 
but it does include non-global logarithms beyond the "one out-of-gap gluon approxi- 
mation" . Another effect which is missing in the parton shower approach is Coulomb 
gluon exchange. As pointed out in [19] these contributions are especially important 
in the large Q/Qo and large Ay region: Coulomb gluons contribute to building up 
the colour-singlet exchange contribution, which eventually leads to a rise of the gap 
fraction at large enough Ay. 

The ATLAS collaboration also compared their data to theoretical predictions 
obtained with HE J [29]. That framework is based on the factorization of multi-gluon 
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amplitudes in the high-energy regime. As in the BFKL approach, a^Ay"^ terms are 
resummed, but energy-momentum conservation is enforced. Logarithms of Q/Qq 
are not systematically resummed, unless they come with a Ay factor. We notice 
that the HE J predictions are similar to ours for the global part, after accounting 
for energy-momentum conservation. This does not come as a surprise: Although 
the two approaches resum different terms, the leading contributions are of the form 
Q/n^yU ij^n ^ ^^^ g^j,g resummed in both approaches. HEJ describes emissions of out- 
of-gap gluons and, if interfaced with a parton shower, should be capable of capturing 
non-global logarithms as well [48]. We note that the HEJ framework does not at 
present include colour mixing via Coulomb gluons. 

It seems clear that within the context of the overall accuracy of a leading 
log/leading order matching and the kinematic range of the current data, the im- 
pact of sub-leading N^, and Coulomb gluon effects is not yet critical (except perhaps 
at the largest values of Ay where both PYTHIA and HERWIG++ undershoot the 
data). The same cannot be said about the constraints of energy-momentum conser- 
vation, which are clearly very important. The message is clear: the accuracy of the 
ATLAS data already demands better theoretical calculations. 
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